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Abstract 



As the ultimate stage of stellar nucleosynthesis, and the source of the iron peak nuclei, 
silicon burning is important to our understanding of the evolution of massive stars and 
supernovae. Our reexamination of silicon burning, using results gleaned from simulation 
work done with a large nuclear network (299 nuclei and more than 3000 reactions) and from 
independent calculations of equilibrium abundance distributions, offers new insights into 
the quasi-equilibrium mechanism and the approach to nuclear statistical equilibrium. We 
find that the degree to which the matter has been neutronized is of great importance, not 
only to the final products but also to the rate of energy generation and the membership 
of the quasi-equilibrium groups. A small increase in the global neutronization results in 
much larger free neutron fluences, increasing the abundances of more neutron-rich nuclei. 
As a result, incomplete silicon burning results in neutron richness among the isotopes of the 
iron peak much larger than the global neutronization would indicate. Finally, we briefiy 
discuss the limitations and pitfalls of models for silicon burning currently employed within 
hydrodynamic models. In a forthcoming paper we will present a new approximation to the 
full nuclear network which preserves the most important features of the large nuclear network 
calculations at a significant improvement in computational speed. Such improved methods 
are ideally suited for hydrodynamic calculations which involve the production of iron peak 
nuclei, where the larger network calculation proves unmanageable. 
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1. Introduction 

With the exhaustion of hydrogen in the core of a massive star, an inexorable contraction 
begins, heating and compressing the core, delayed for a time as each succeeding nuclear fuel 
ignites and is transformed. Beginning with the helium ash of hydrogen burning, most of these 
burning stages consist of fusion reactions among the nuclei of the ash. The first exception 
is neon burning. Following carbon burning, which leaves behind ash composed of O, Ne, 
and Mg, the temperature and density continue to rise. Before the temperature is adequate 
to allow fusion reactions among O nuclei to overcome the Coulomb repulsion, the photon 
field becomes sufficiently energetic to photodissociate Ne. In general, a photodisintegration 
channel becomes important when the Q-value of the reaction, that is the energy difference 
between fuel and products, is smaller than approximately SO/c^T. For such temperatures, the 
high energy tail of the Planck distribution provides enough photons of sufficient energy that 
the photodisintegration reaction represents a comparable flow. With a Q-value of 4.7 MeV, 
^°Ne(7, a)^^0 becomes important for temperatures above 1.5 x 10^ K. Thus the next stage is 
neon burning, typified by the photodisintegration of neon and the subsequent capture of the 
ejected a-particles by the remaining heavy ions. With the exhaustion of this energy source, 
the collapse continues until the O nuclei are sufficiently energetic to fuse. Once more nuclear 
burning offers a respite until the oxygen fuel is exhausted and the collapse continues, further 
raising the temperature and density. With typical Q- values for reactions among stable nuclei 
above silicon being 8-12 MeV, photodisintegration begins to play an important role for all 
nuclei once the temperature exceeds 3 x 10^ K. Nuclei with the smallest binding energies arc 
destroyed in favor of their more tightly bound neighbors. Very quickly this leaves the core 
composed of silicon isotopes and their a-nuclei neighbors, S, Ar, and Ca, the most bound of 
the light nuclei. Continued contraction further increases temperatures, favoring even more 
tightly bound nuclei, the iron peak. What follows is a complex scries of photodissociation 
and capture reactions, converting silicon into iron peak elements. Unlike previous burning 
stages, where a few reactions (either heavy ion reactions among the principal constituents 
or conversion among neighbors initiated by photodisintegration) dominated, here the fuel 
nuclei are linked to the product nuclei by a multitude of reaction chains and cycles. Since 
these chains of reactions wind their way through many nuclei intermediate between silicon 
and the iron peak nuclei, it is necessary to keep track of many more nuclei than was the 
case for prior burning stages. With four potential particles in the incoming channel, p, n, 
a, or 7, and these same four choices for the outgoing channels, discounting elastic scattering 
reactions, there are, in principle, at least 12 reactions per nucleus. Thus silicon burning is a 
complicated web of reactions, making it necessary to keep track of the abundances of a large 
number of nuclei. This is a major contribution to the complexity of silicon burning. 

Silicon burning is further complicated by the closeness to equilibrium of many of these 
pairs of photo dissociations and capture reactions. Net fluxes are often orders of magnitude 
smaller than either the forward or reverse reaction rates would indicate. Considering that 
the end state of silicon burning is an equilibrium, nuclear statistical equilibrium (NSE), it is 
hardly surprising that equilibria among nuclei arise during silicon burning. From a purely 
numerical point of view, the net reaction flux being the difference of two large numbers 
raises the potential danger of round-off errors. However, equilibria arc attractive for the 
simplification they provide. As we will see, equilibria are the key to understanding silicon 
burning and simplifying our modeling thereof. 

The work of Bodansky, Clayton & Fowler (1968; hereafter BCF) examined the process 
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by which this equihbrium distribution is formed. Burbidge, Burbidge, Fowler & Hoyle (1957; 
hereafter B^FH) had postulated that the intermediate nuclei, principally the a-particle nuclei 
from 24Mg to ^^Ca, were formed by the a-process, a series of a particle captures, which was 
separate from the e-process which formed the iron peak. Taking their cue from Hayashi et 
al. (1959) and the nuclear network calculations of Truran, Cameron, & Gilbert (1965), BCF 
showed that the intermediate nuclei are formed via a partial equilibrium in which groups 
of nuclei exist which are internally in equilibrium under the exchange of photons, protons, 
neutrons and a-particles, i.e. in equilibrium with respect to strong and electromagnetic 
reactions. These quasi-equilibrium (QSE) groups remain out of equilibrium with respect 
to other groups until the final complete equilibrium is achieved. BCF concluded that the 
abundances of the intermediate nuclei could be explained by a single QSE group which 
reached from ^^Si through the iron peak, but which failed to reach complete equilibrium, or to 
completely exhaust sihcon. Thus the a- and e processes of B^FH were unified, coinciding with 
silicon burning, which had been determined to be a fundamental late stage of nuclear burning 
in massive stars. Further, BCF examined the influence of weak reactions on silicon burning 
and concluded that, while these reactions could have important effects on the abundance 
distribution, the ratio of the total number of protons to the total number of neutrons did 
not vary far from unity. In contrast to Fowler & Hoyle (1964), BCF concluded that the 
dominance of ^^Fe was due to the production of ^^Ni which subsequently decayed to ^^Co 
and then to ^^Fe. This scenario agrees well with observations of supernovae. The lack of an 
equilibrium among weak reactions requires the monitoring of an additional degree of freedom 
which affects the equilibrium solution, the degree to which the material is neutronized. In 
the literature this is parameterized in two ways, as Yg, the electron (molar) abundance, or 
as 77, the neutron excess parameter: 

i i ^ *^ 

,J = J2[Ni-Zi]Y^ = J2 

i i 

V^i-'2Ye, (Ic) 

where Zj, iVj, Ai, Xi, Yi are the atomic number, neutron number, mass number, mass frac- 
tion, and abundance of nucleus i. Physically, Yg is the ratio of protons to nucleons (identical 
to the ratio of electrons to nucleons, hence the subscript e) or the total proton fraction, and 
rj is the fraction of excess neutrons per nucleon. BCF showed that the abundances which 
result from silicon burning are strongly dependent on the degree of neutronization, a view 
reinforced by Hartmann, Woosley, & El Eid (1985). In this paper we will show that not 
only the nuclear products but the entire mechanism is strongly affected by the degree of 
neutronization. 

Woosley, Arnett & Clayton (1973; hereafter WAC), in the context of parameterized ex- 
plosive silicon burning, showed that rather than a single QSE group between silicon and 
the iron peak, there are initially two groups, roughly divided by A ~ 45. WAC further 
demonstrated that considering only the reactions which link the groups yields a good ap- 
proximation of the time necessary for the two groups to merge and form a single QSE group. 
Indeed, WAC contended that this linkage was dominated by a single reaction, ^^Sc(p, 7)^^Ti, 
with perhaps a quarter of the flow going through less important reactions like ^^Ca(Q;, 7)^^Ti 



Ni-Zi 



Ai - 2Zi 



Xi 



(lb) 
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and Ti(n, 7) Ti. In later sections, we will show that for material which has been more 
highly neutronized, these two QSE groups form further from mutual equilibrium and that 
this division persists for a much longer time. In addition there is some ambiguity in defining 
the boundary between the QSE groups and the most important reactions linking the groups 
as Ye is varied. As we will also show, the dependence of the behavior of the QSE groups on 
the degree of neutronization has important consequences for the energy generation and other 
physical manifestations of silicon burning, in addition to determining the nuclear products. 

Thielemann & Arnett (1985; hereafter TA), examined silicon burning in the context 
of hydrostatic models of massive stars and noticed behavior largely in keeping with that 
described by WAC for the explosive case. For conditions characteristic of the cores of more 
massive stars, high temperature, low density, and consequently larger Ye, TA found that the 
bottleneck between the QSE groups, coinciding roughly with Z = 21, was bridged on the 
proton-rich side of stability. However, for conditions more characteristic of lower mass stars, 
lower temperature, higher density, and smaller Ye, the bridge was found to be proton capture 
on neutron-rich isotopes of Ca. This more complex behavior differed with the assertions of 
WAC and of Weaver, Woosley, & Fuller (1985), that the single reaction ^^Sc(p, 7)^^Ti is 
responsible for the lion's share of the flow between the QSE groups. As we will show in later 
sections, our efforts allow us to explain the differences between the results of TA and WAC. 

2. Nuclear Reaction Network Calculation 

Silicon burning in the cores of massive stars takes place under a range of temperatures and 
densities and with a range of electron fractions. We have examined silicon burning in a cube 
of this parameter space with temperatures in the range 3.5 to 5 xlO'^'K, densities between 
10^ and 10^° gem , and Yg between .498 and .46, using a nuclear network of 299 nuclei, 
shown in Tabic 1. This region of parameter space roughly spans the region that previous 
investigations of nucleosynthesis in massive stars (Thielemann & Arnett 1985, Nomoto & 
Hashimoto 1988) have exhibited during core silicon burning. Further, the lower density 
portions of this parameter space overlap, to a large extent, the parameter space that prior 
investigations (WAC, Woosley, Pinto, & Weaver 1988, Thielemann, Hashimoto, & Nomoto 
1990, Aufderheide, Baron, & Thielemann 1991) have determined appropriate for explosive 
silicon burning. We have included the effects of Coulomb screening on the equilibria for 
this same region of parameter space, as discussed in Hix, Thielemann, Fushiki, & Truran 
(1996). The results of nuclear network calculations for this region of parameter space and 
the comparison of these network calculations to equilibrium calculations are discussed in 
the present paper. In a subsequent paper (Hix & Thielemann 1996, hereafter Paper II), 
we will discuss the applicability of quasi-equilibrium to explosive burning. The final goal is 
an improvement to the network, using what we have learned about quasi-equilibrium, which 
speeds the calculation of the energy generation and nuclear abundance changes due to silicon 
burning, preserving much of the accuracy of the network calculation while greatly reducing 
the computational overhead. 

For our study of silicon burning, we consider 299 nuclei, listed in Table 1, linked by more 
than 3000 reactions. This nuclear set stretches from protons and neutrons to germanium. 
With 11 or more isotopes per element around and above iron and 7 or more in the region 
around silicon, this nuclear set is more complete than that used in earlier work discussed in 
§1. This improvement is most important for material which has undergone significant prior 
electron capture and hence has a Ye significantly less than .498. Rather than complicate 
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Table 1: Nuclear set Included in Calculations 
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our task by including hydrodynamics, we have instead done a parameter study of silicon 
burning as a function of three variables, temperature, density and Y^.. In the course of 
a hydrodynamic model calculation, the star's core would wind its way through our cube 
of parameter space. In particular, electron capture will cause Kg to decrease with time. 
However, the timescales for changes in Fg via weak reactions, even at densities sufficient 
for significant electron capture, are much longer than the timescales for the strong and 
electromagnetic reactions. Thus it is possible to treat them separately. Therefore we have 
neglected weak reactions, allowing us to treat Ye as a constant parameter rather than merely 
considering the Ye at one point in time. The reaction rates for light nuclei are based on 
experimental information from Caughlan & Fowler (1988), Bao & Kappeler (1987), Wiescher 
et al. (1986, 1987, 1989), Wiescher, Gorres, & Thielemann (1988), and also from Wagoner 
(1969) and Wagoner, Fowler, & Hoyle (1967) if otherwise unavailable. For intermediate and 
heavy nuclei, where higher level densities permit application of statistical model calculations, 
rates by Thielemann, Arnould, & Truran (1987) are used (see also Cowan, Thielemann, & 
Truran 1991). We utilized a reaction network at constant temperature, density, and Yg, 
with Ye set by an initial distribution of silicon isotopes. There are of course several such 
distributions with the same Yg but, with the exception of a brief initial adjustment phase 
which does not enter into our analysis, such variations of the initial distribution showed the 
same results. The basics of a nuclear reaction network are summarized in Woosley (1986) 
or Thielemann, Nomoto, & Hashimoto (1994). It is particularly important to note that the 
weak and intermediate screening prescriptions of Graboske et al. (1973) were used, as well as 
the strong screening prescription of Itoh ct al. (1990). Transition among these prescriptions 
is performed by utilizing the smallest predicted screening enhancement. 

3. The Physics of Quasi-Equihbrium 

Previous authors, notably BCF and WAG, have shown that quasi-equilibrium (QSE) is 
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a most important aid to understanding the process of silicon burning. This is fortunate as it 
is much easier to follow the evolution of equilibrium groups rather than become lost in the 
welter of individual abundances and reactions. Since ^^Si is the principal fuel, we will begin 
our discussion of QSE here. If we take ^^Si to be the focus of our quasi-equihbrium group, 
then the abundance of a nucleus which is in equilibrium with ^*Si with respect to the 
exchange of free nucleons, ct-particles, and photons is 



VC(28Si) 



('i&Q^Wv N-\^^^ Z-U 



(2) 



where F^, Yp, and F(2^Si) are the abundances of free neutrons, free protons, and '^^Si, and 



28; 



we have defined 



B{^Z) 



(3) 



for later convenience. G{^Z) and B{^Z) are the partition function and binding energy of 
the nucleus ^Z, Na is Avogadro's number, ks is Boltzmann's constant, and p and T are 
the density and temperature of the plasma. Thus, the abundance of a nucleus in quasi- 
equilibrium with ^^Si is a function of three other abundances, those of free protons, free 
neutrons, and ^^Si, and the thermodynamic conditions and properties of the nucleus, sub- 
sumed here within C{^Z). It is the evolution of these three abundances which determine 
the behavior of the entire QSE group. We reserve Yqse{^Z) to represent the abundance of 
a nuclear species in quasi-equilibrium with ^^Si. Equation (2) is identical to the expressions 
for quasi-equilibrium introduced by BCF and WAC, provided the a-particles, protons and 
neutrons are internally in equilibrium. We find this provision to be justified by the time 
quasi-equilibrium is established. 

This derivation ignores screening, an omission which Hix et al. (1996) found can be 
significant for equilibria under the conditions we are discussing. A rederivation including the 
screening corrections to the reaction balance produces 



VC(28Si) 



Z-1 

n 

Z'=U 



exp{Hp,y{Z')) 



(4) 



where cxp{Hij{Z')) is the screening correction for the reaction Z'{i,j)^ Z" , and we utilize 
the convention that, for indices less than the bottom index [Z < 14), the product is instead 
a division. The form of this equation is perhaps not surprising. We can rewrite Equation 
(4) in the suggestive form 



C{^Z)Yn^Yp^Ul-2i 



exp(gp,^(ZO) - 

exp(ff^,p(Z')). 



F(28Si) 



C{^^Si)Yn'%'^Uz'=i 



(5) 



As F(28Si), Yp, and Y^ approach their NSE values. 
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r(28si) ^ Y^sEe'si) ^ ce'si)Yr,'%'^ n 



Z'=l 



' exp{Hp^^{Z')) 



(6) 
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forcing 



Yqse{^Z) C{^Z)Yn''Y/ 



Z-l 

n 

Z'=l 



(7) 



Thus the quasi-equihbrium group blends simply into the NSE distribution as 28Si and the 
free nucleons approach equilibrium. 

To compare the network results to quasi-equilibrium, we define 



rQSE{^Z) = log 



QSE[ 



(A 



Z) 



(8) 



with Y(-^Z) being the network abundance; Yqse{^Z) is calculated from the network abun- 
dances of free nucleons and ^^Si. For comparison, this definition is compatible with Vqe as 
defined by WAC. Figures la-g show the results of this comparison for varying degrees of 
silicon exhaustion, at different temperatures, densities, and values of Y^. Before we consider 
the results, it is instructive to determine the signature that multiple quasi-equilibrium groups 
leave on rgsEi^Z). 

If a nucleus "^Z is in quasi-equilibrium with, for example, ^^Ni, then there is an expression 
analogous to Equation (4) or Equation (5) for this abundance. Taking Equation (5) as the 
template, we can define a second quasi-equilibrium abundance with respect to ^^Ni, 



Y 



QSE2[ 



C{^Z)Yn''Y/Ylz'A 



Yi^^Ni) 



c{^^m)Y,,'%'^U'z'=i 



(9) 



For species which are part of this second QSE group, i.e., Y(/^Z) — Yqse2{^Z), using 
Equation (9) and Equation (5) we can evaluate 

Yqse{^Z) Yqse{^Z) 



Ny ZyfZ-l \ cy.p{H.,,~,(z')) ' 
I' iiz'-i ,.._xiHJi-,.,,(Z')) 



C7(28si)y„i4yp"n?'=i 



cxp(gp,^(ZO) 
exp(i/^,p(Z')) 



Y{^Z) YqsE2{^Z) f CiAZ)Y^VUz", 



exp(gp,^(ZO) 

oxp(H-/,p(Z')) 



(10) 



C(56Ni)y„28y^28-Q^7^^ 

Since the numerators of these fractions are identical, cancellation leaves 



exp(gp,7(Z0) 
exp(H^,p(Z')) 



r(56Ni) 



YqseC'Z) I 


{ce^m)Y,,^%^'UlU 


rexp(ii;,.-,(Z'))" 
[cxp(i/^,p(Z'))_ 
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Thus every member of a quasi-equilibrium group has the same value of rqsEi^Z). For the 
silicon group, rQSE{^Z) = 0, but other groups will also be apparent but offset from zero. 
Although free nucleons do not come into equilibrium with ^^Si until NSE is established, it 
is also instructive to examine rqsEip) or rqsEip)- From Equations (3) and (5) we have 
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Figure la: Comparison of the silicon quasi-equilibrium abundance to the network abundance, as a function 
of A for Tg = 5.0, p = lO^gcm-^, and Yg = -498, with X(Si group)=.99. 



Dividing by Yp gives 



y(28si) 



P C{^^Si)Yn'%''nz'=i 



Yr,. 



(13) 



The denominator of Equation (13) is identical to the expression for the NSE abundance of 
^^Si, reflecting the abundance of ^^Si if it was in equilibrium with the free nucleons. Thus 
changes in rQSE{p){= ^'QSE{n)) reflect the approach of the silicon group to NSE. 

4. Compcirison of Network results to Quasi-Equilibrium 

In Figures la-g we compare the results of the full network with abundance predictions 

calculated from the network abundances of n, p, and ^®Si assuming quasi-equilibrium. El- 
ements in quasi-equilibrium will show similar values of rq^^i^Z) in each of Figures la-g. 
Figure la shows what is a representative pattern when rqsEi^Z) is plotted against A. From 
roughly A=24 to A=45, there is a large cluster of nuclei in quasi-equilibrium with ^^Si, with 
a second quasi-equilibrium group stretching from roughly A=50 to the top of the network. 
Finally, there is a third small group composed of protons, neutrons, a-particles and some 
other light nuclei. Species in between this lightest group and the silicon QSE group do not 
form large QSE groups, although they are closer to quasi-equilibrium with silicon than the 
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Figure lb: Comparison of the silicon quasi-equilibrium abundance to the network abundance, as a function 
of A for Tg = 5.0, p = lO^gcm-^, and Yg = -498, with X(Si group)=.90. 



light element group. With a temperature of 5 x 10 K, a density of 10 gcm~ and Ye of 
.498, this case is among the fastest we will consider. After an elapsed time of 3.0 x 10~^ s, 
with 86% of the mass contained in ^^Si (X(28Si)=.86), and 99% of the mass in the silicon 
QSE group (X(Si group)=.99), quasi-equilibrium is clearly well established. Comparison 
with Figure lb shows that at a later time (elapsed time = 5.6 x 10~^ s, X(Si group)=.90, 
X(^^Si)=.65), this pattern is maintained although all nuclei are closer to equilibrium. In fact 
all of the nuclei above Na are within 25% of their sihcon QSE abundance. However exami- 
nation of Figure Ic, a detail of Figure lb, plotted now as a function of Z, shows that even 
within this narrow margin it is still possible to resolve the 2 QSE groups and a scattering 
of nuclei between them. Since rqsEi^Z) is plotted against Z in this figure, with important 
nuclei labeled by their respective atomic masses, it is also immediately apparent that the 
boundary between the groups is not as simple as Z = 21, since neutron- rich isotopes of K 
and Ca are clearly members of the upper QSE group while proton-rich isotopes of V and 
Cr are in quasi-equilibrium with Si. The description of the boundary as A = 45 is more 
successful, giving way to N=23 from Ti onward. Examination of Table 2 allows a more quan- 
titative comparison of the values of rQSEi^Z), a point we will return to in a moment. The 
convergence toward equilibrium continues as time elapses. By an elapsed time of 2 x 10^^ s 
(X(Si group)=.5), the elements that formed the upper group are within 5% of their silicon 
QSE abundance, even nuclei with mass fractions as small as 10~^^. 



10 



6 

o 



.15 



.05 



-.05 



-.1 



T9=5.0 
Ye=.498 



p=iO g cm 
Time=5.61xl0"^ s 



V Cr Mn Fe Co Ni 



46° 11° il^ 468 478 48° 49o 50' 
450 470 48o 

44° 440 440 

I I 438 438 "^"^0 441 458 46° 
Ar K Ca Sc Ti 



S i -i 



10 12 14 16 18 20 22 24 

Z (Atomic Number) 



26 



28 



30 



32 



Figure Ic: Detail of the comparison of the siHcon quasi-equihbrium abundance to the network abundance, 
as a function of Z for Tg = 5.0, p = 10^ gcm"^, and Yg = -498, with X(Si group)=.90. 



Previous authors have discussed the progress of sihcon burning in terms of the degree 
of exhaustion of ^^Si. For this investigation, where Ye is a parameter, this description is 
suspect. The abundance of ^®Si in a composition purely of sihcon decreases as Yf. decreases. 
In place of X(^®Si), we have chosen to use the mass fraction within the entire silicon QSE 
group, X(Si group). This choice has several advantages. Since the mass fraction is dominated 
by the two QSE groups, the mass fraction within the sihcon group does reflect the degree to 
which material remains to be passed into the upper group (which in NSE will compose 95% 
or more of the mass). Also, this choice is useful for all in a way not possible for a single 
nucleus. This choice is, however, limited in the sense that individual reactions depend on 
the abundances of individual nuclei. Knowledge of the mass fraction in the silicon group is 
not sufficient to determine the photodisintegration rate of ^^Si, for example. This is an issue 
which we must keep in mind as we proceed. 

As a further demonstration of this difficulty, consider the case where Tg = 5.0, p = 
lO^gcm"^, but Ye = .46. Figure Id is comparable to Figure lb, in the sense that X(Si 
group) is approximately the same (.91 for Figure Id with Yg = .46 as compared to .90). 
But the fraction of ^^Si differs greatly, with X(^^Si)=.008 for Ye = .46. Furthermore, for 
Ye = .46, X(28Si) is actually rising slightly at X(Si group) ~ .9. It should be noted, however, 
that for Ye = .46, X(all isotopes of Si) is .64, thus this extremely small abundance of ^®Si is 
specific to this nucleus, and is due to the preference for more neutron- rich nuclei at lower Ye- 
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Figure Id: 
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Comparison of the silicon quasi-equilibrium abundance to the network abundance, as a function 
= 5.0, p = lO'^gcm-3, and = .46, with X(Si group)=.90. 



There is also considerable difference in the elapsed time, 1.1 x lO^"' s for Yg = .46 compared 
to 5.6 X 10~^ s for Yg — .498. Although these elapsed times are small, they are significant, 
representing the time required to consume 10% of the silicon fuel. Thus such differences 
in elapsed time are important, implying faster transfer of mass to the iron peak group and 
consequently, as we will discuss in §7, a larger energy generation rate. For now however, the 
concern is quasi-equilibrium. Comparison of Figure Id with Figure lb shows that the gap 
between the silicon group and the iron peak group is much larger for Yg — .46. Indeed, the 
underabundance of the iron peak group is more in keeping with that shown in Figure la. 
Comparison of Figures Ic and Ic reveals that, rather than being 25% undcrabundant, the 
nuclei in the QSE group around iron are underabundant by a factor of 25 for Y^ = .46. This 
is in spite of the similar fractions of the mass which have been transferred to the iron peak 
group, .09 and .1 for Y^ — .498 and .46, respectively. Clearly, the relative underabundance of 
the iron peak group at low Ye reflects a much larger QSE abundance for these species rather 
than an actual dearth of mass transferred. Another point to note in Figure Id is the larger 
value of TQSEip) for Yg — .46 compared to Yg — .498. This reflects the smaller NSE mass 
fraction of the silicon group (.00005 for Ye = .46, compared to .018 for Yg = .498), for lower 
Ye and hence for the same degree of fuel exhaustion, the composition is further from NSE. 
These arguments indicate pointedly that the quasi-equilibrium behavior of silicon burning 
is strongly influenced by the state of neutronization of the material. 
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Figure le: Detail of the comparison of the silicon quasi-equilibrium abundance to the network abundance, 
as a function of Z for Tg = 5.0, p = 10''gcm~'^, and Yg = .46, with X(Si group)=.90. 



Another illustration of the impact of neutronization, which is not easily discernible from 
Figures lb and d, but is very apparent in Figures Ic and e, is the changing constituents 
of the QSE groups with declining Yg. For Yg = .498 the distribution is well described as 
2 groups, separated by a narrow fringe of nuclei with A = 45 up to Ti and then = 23. 
The similarity of the value of rqsEi/^Z) for those 4 species with A" = 23 (^^Ti, ^^V, ^^Cr 
and ^^Mn) implies that the reactions hnking them are close to being balanced. Comparison 
of Figure Ic and le reveals a much more complicated pattern. In addition to the vastly 
different scale, there are many more nuclei in the boundary region. Instead of only a pair of 
Ca isotopes (^"^Ca and ^^Ca), the isotopes from '^^Ca to ^^Ca are all significantly displaced 
from the QSE groups. A similar situation exists for Sc and Ti. Those isotopes of Ti, V, 
Cr, and Mn which were in the silicon group or the boundary for Y^ — .498 are significantly 
closer to the upper group. The nuclei with the largest A which are within a factor of 2 of 
quasi-equilibrium with the silicon group are ^2Sc and ^^K. The A^ = 23 nuclei are much 
closer to the upper group, with ^^Cr, for example, having an abundance 20% larger than its 
iron peak group QSE value as opposed to being a factor of 20 underabundant for QSE with 
the silicon group. From Ca on, the N — 22 and A^ = 21 elements occupy the boundary. The 
vertical structure of these lines of constant A^ implies that the reactions linking them, pairs 
of (p, 7), (7,p) reactions, are not in equilibrium. Similarly the spacing of isotopes of Ca, for 
example, indicates an unbalanced fiow by neutron capture upward from the silicon group. 
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We will examine in detail the reaction flows between the groups in §6. One more interesting 
note is that, in contrast to the behavior of other nuclei, the neutron-rich isotopes of K (^'^K, 
^^K, and ^^K), although actually further from their silicon group QSE abundances, seem to 
be closer to membership in the silicon group. Table 2 shows quantitatively the movement of 
nuclei upward from the silicon group, with those abundances below the minimum resolution 
of the network (10~^^) represented by •. Comparison of the second and third columns shows 
clearly the migration of the boundary toward lower Z and A. The abundance of a nucleus 
like ^^Sc, which for = .498, is within .5% of its QSE abundance and thus clearly a member 



14 



-.1 



T9=5.0 

Ye=.46 



p=10' g cm 
Time=5.60xl0 s 



48-490 



V Cr Mn Fe Co Ni 



% 46o 



s e 4/9 



46-478 ^§8 



43o 44o 



45-46§ 
43-44g 



420 



430 



46o 



450 



440 



46o 



I e i 



440 



42o 



42o 



e i i 1 1^^1418 418 

Ar K Ca Sc Ti 



-.2 



10 12 14 16 18 20 22 24 

Z (Atomic Number) 



26 



28 



30 



32 



Figure If: Detail of the comparison of the silicon quasi-equilibrium abundance to the network abundance, 
as a function of Z for Tg = 5.0, p = 10^ gcm~^, and Yg = .46, with X(Si group)=.52. 



of the silicon group, is, at Yg — .46, better predicted by QSE with the iron peak group. For 
Ye = .46 and X(Si group) ~ .9, the sihcon QSE abundance of ^^Sc is over 11 times its network 
abundance, while its iron peak group QSE abundance is slightly less than half its network 
abundance. And ^^Sc is clearly not an isolated example. Thus the group boundaries defined 
by earlier work were clearly dependent on the narrow range of parameters examined. We 
will therefore investigate how our broader range in parameter space affects these boundaries. 

One possibility is that the differences in the boundaries are due to the fact that the 
lower Ye cases are further from equilibrium. This would result in a temporal variation in the 
boundary with the boundary discussed above for lower Ye approaching that of the Ye = .498 
case at X(Si group) ~ .9, as the lower Yq cases reach a similar degree of convergence. We 
find that this is not the case. Consider the situations for X(Si group) ~ .5. This corresponds 
to elapsed times of 2.1 x 10~^ s and 5.6 x 10~^ s, for Ye = .498 and .46, respectively, with 
Tg = 5.0 and p = lO^gcm""^. As the fourth column of Table 2 demonstrates, for Ye = .498 
the abundances of all of the species which composed the iron peak QSE group are within 
2-3% of their sihcon quasi-equihbrium abundances. Figure If illustrates the comparison of 
the network abundances with those predicted by quasi-equilibrium for Yg = .46 and X(Si 
group) ~ .5. The members of the iron peak group are within 50% of their silicon quasi- 
equilibrium abundances. Comparison of Figures Ic, le, and If reveals nonetheless that the 
pattern for Y^ — .46, X(Si group) ~ .5 is more like that of Y^ — .46, X(Si group) ~ .9, 
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Figure Ig: Detail of the comparison of the silicon quasi-equilibrium abundance to the network abundance, 
as a function of Z for Tg = 3.5, p = 10^ gcm'^, and = .498, with X(Si group)=.90. 



even although the iron peak group has converged much closer to quasi-equilibrium with Si. 
Although the increased resolution reveals much more detail among the upper members of the 
boundary, and the isotopes of K continue their slow convergence toward the silicon group, 
a greater degree of silicon exhaustion leaves the pattern essentially unchanged. While the 
break which leaves '^^Sc and ^^Ti, newly emerged from below the minimum mass threshold 
of the network, close to the silicon group is made more prominent by the approach of these 
nuclei to the silicon group, the constituents of the QSE groups are essentially unchanged. 
Furthermore, the dissimilarity of Figures Ic and If, both corresponding to elapsed times of 
5.6 X 10^^ s, removes any doubt that the differences in the group boundaries are simply a 
time dependency. Thus we can conclude that the differences in the boundary between the 
QSE groups is a genuine effect of the change in Yg and not due to temporal variation or 
closeness to equilibrium. As we will discuss in §7, with the Yg = -46 case having exhausted 
5 times as much fuel as the similar case with Fg = -498 in the same elapsed time, these 
differences in QSE group structure have broader implications than simply the abundances 
of these boundary nuclei, resulting in large variations in the rate of energy generation. 

Having shown that the quasi-equilibrium behavior during silicon burning is dependent 

on Ye and the degree of silicon exhaustion, we now turn to our two remaining parameters. 
The conditions considered so far, high temperature and low density, are the most like those 
considered by WAC. They are also the conditions under which the inclusion of screening is 
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the least critical. For conditions like those of Figure Ig, Tg = 3.5 and p = lO^gcm"'^, the 
inclusion of screening in the quasi-equilibrium is very important. Under these conditions the 
screening portion of Equation (10) can be as large as 10^, much larger than the factor of 5 
which is the maximum for Tg — 5.0 and p — lO^gcm"^, but much smaller than the factor 
of order 10^^ which is possible for Tg = 3.5 and p — lO^^gcm"^. Comparison of Figure Ic 
and Ig supports the inclusion of this screening term. With Ye = .498 and X(Si group) = 
.9, the conditions of Figure Ig differ from the conditions of Figure Ic only in temperature 
and density. Since the screening term is dependent on Z and not on N, its omission would 
result in a separate horizontal line for each element, with Ge being displaced from Si by a 
factor of almost 10^. The similarity of Figure Ic and Ig argues strongly that the inclusion 
of screening is correct. This similarity also argues that changes of temperature and density 
do not strongly affect the structure of the QSE groups, although the elapsed time. 111 s, is 
considerably different, reflecting the much slower reaction rates at low temperatures. The 
most striking difference between these cases is the larger values of VQSEi^Z) among the 
lighter nuclei. The much greater value of rgssip) (16.92 compared with 7.32) for this low 
temperature/high density case is consistent with the argument that rqsEip) reflects how far 
from NSE the distribution is. Since the NSE mass fraction of the sihcon group for these 
conditions is less than 2 x 10~^, the smallest of the cases shown in this paper, one would 
expect rgssip) for this degree of silicon exhaustion to be the largest in this case. 

Clearly changes in temperature and density affect the reaction rates and hence timescales, 
but how does this affect adherence to the quasi-equilibrium distributions? Comparison of 
Figure Ig with Figure Ic shows that, for Yg = .498, the variation of temperature and density 
has little effect on the QSE group behavior. Once again, for X(Si group) = .9, the members 
of the iron peak group are within 25% of their silicon QSE group abundances. The spread 
within the iron peak group is noticeably larger, but still less than approximately 5%. Further 
investigation reveals that this greater disorder, i.e. larger variation within the QSE groups, 
is solely a function of decreasing temperature. Most importantly, although there is some 
difference in relative placement, the nuclei intermediate between the two QSE groups are 
the same. Examination of cases with lower Yg supports the finding that the variation of 
temperature and density does not afi^cct the grouping of nuclei, although it can radically 
alter relative abundances among the collection intermediate to the QSE groups. 

While previous authors have shown the value of quasi-equilibrium as a key to under- 
standing the process of silicon burning, the work presented in this section implies that this 
understanding is incomplete. As one might expect, variation of the temperature or den- 
sity can produce large differences in the behavior of silicon burning, most notably in the 
timcscale. The important result of this research is the importance of neutronization in the 
quasi-equilibrium behavior and timescale. For similar degrees of silicon exhaustion, a de- 
crease in Ye results in a drastic increase in the QSE abundances of the iron peak group. 
As a result, for similar degrees of silicon exhaustion the iron peak group is much farther 
from quasi-equilibrium and the abundances of these nuclei a much smaller fraction of their 
silicon QSE abundance. This delays the merging of the two QSE groups into a single group 
spanning the elements from Mg to Ge. While high Ye cases reach this merged group stage 
with approximately 90% of the mass remaining within what was the silicon group, for lower 
Yf, this mass fraction is much smaller. For Y^ = .498, for both Tg = 5.0, p = lO'^gcm^'^ and 
Tg = 3.5, p = lO^gcm""^, a single QSE group (within 10%) has formed by X(Si group) = 
.85. For Ye = -48 the iron peak nuclei do not reach 90% of their silicon QSE group abundance 
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until 50% or more of the mass in the sihcon group is exhausted. At Yg = .46 for Tg = 5.0, 
p = lO^gcm"^ and Tg = 3.5, p = lO^gcm^"^ the iron peak group reaches 90% of its sih- 
con QSE abundance with sihcon group mass fractions of .25 and .3, respectively. Although 
variations of temperature and density clearly effect great changes in the timescales of silicon 
burning, it is the degree of neutronization which most strongly affects the quasi-equilibrium 
behavior. This behavior, previously only hinted at by TA, is very important to modeling 
the mechanism of silicon burning, since low Yg results in the persistence of 2 separate QSE 
groups until much later than previously thought. Attempts, like that of BCF, to model 
silicon burning by using the quasi-equilibrium abundances and calculating the loss from a 
single QSE group are, as a result, oversimplified. Our results show that such a method would 
be inaccurate unless this separation between the upper and lower QSE groups is included, 
as the gap can persist through significant silicon consumption, causing large miscalculations 
in the iron peak abundances. 

Furthermore, for smaller Yg, the upper boundary of the silicon group moves downward, 
while at the same time the nuclei which formerly comprised the boundary layer have, for 
the most part, joined the iron peak group. This movement of the boundary has a number 
of effects, the most important being that it mandates there be different reactions which link 
the QSE groups as Y^ varies. It is hard to see how the reaction ^^Sc(p, 7)^^Ti can be the 
link between QSE groups for Yg = .46, since ^^Sc is virtually a member of the iron peak 
group under these conditions. The similarity in r QSEi^Z) for ^^Sc and ^^Ti for Yg = .46 
implies that the (p, 7), (7,p) reaction pair linking them are almost balanced, making this an 
unlikely flow path. In §6 we will investigate the question of how the linking reactions vary 
as conditions change, but first we will use the convenience of quasi-equilibrium to follow the 
nuclear abundances produced during silicon burning. 

5. Species Formation during Silicon Burning 

The applicability of quasi-equilibrium greatly simplifies the study of the nuclei produced 
by silicon burning. Instead of following the intricacies of the production of 300 nuclei, much 
can be learned by examining four: free protons, free neutrons, and one member of each of 
the QSE groups. While the focal members of the QSE groups map the spacing between the 
groups, it is the abundances of the free particles which control the relative abundances within 
the QSE groups. The gap between QSE groups, seen in Figures la-g, which diminishes as 
silicon approaches exhaustion, reflects the evolution of the focal elements of the QSE groups. 
Figure 2 shows samples of the evolution of the free proton and free neutron mass fractions 
as a function of time. Although the actual abundances vary by many orders of magnitude, 
comparison with the appropriate NSE abundances reveals a similarity in pattern. For each 
of the cases portrayed in Figure 2, Tg = 5.0, p = lO^gcm^"^, and Yg = .498 (solid lines), 
Tg = 5.0, p = lO^gcm-^, and Yg = .48 (dotted lines), and Tg = 3.5, p = lO^gcm-^, 
and Yg = .46 (dashed lines), the pattern is similar. Although the rates of convergence vary 
greatly, in each case the network abundance of free protons is initially more than an order of 
magnitude smaller than its equilibrium abundance and converges gradually to equilibrium 
as the mass fraction of the silicon group approaches its equilibrium value. The behavior of 
the neutrons is slightly more complicated. These abundances begin more than an order of 
magnitude larger than their equilibrium abundance, and converge down toward equilibrium, 
briefly becoming underabundant before reaching equilibrium. As with the protons, the 
details of this behavior; the rate of convergence, the maximal and minimal abundances, and 
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Figure 2: Convergence of the free nucleon abundances to equilibrium as a function of degree of silicon 
exhaustion. The cases portrayed are Tg = 5.0, p = lO^gcm"^, and = .498 (solid lines), Tg = 5.0, 
p = 10^gcm~3, and Ye = -48 (dotted lines), and Tg = 3.5, p = lO^gcm"^^ ^nd Yg = -46 (dashed lines). 

the equilibrium mass fractions, vary widely. But this simple pattern is maintained. What 
effect do such variations have? The relative abundance of two isotopes of a given element 
are related by 

Thus an overly large neutron abundance tips the the production of an element toward more 
neutron-rich isotopes. Clearly this occurs because the increased neutron abundance increases 
the rate of neutron capture. Similarly the dearth of free protons will favor the enhancement 
of the members of a QSE group with lower atomic number. Figure 3a shows the evolution 
of the abundance for key members of the iron peak group as the degree of silicon exhaustion 
mcreases for Tg = 5.0, p = lO^gcm'3^ ^nd Ye = .498. For comparison, the corresponding 
evolution of the free nucleons is shown as the solid line in Figure 2. At early time, (X(Si 
group ~ .9)), the iron peak group is dominated by ^^Fe and ^^Fe. As X(Si group) decreases, 
the abundance of ^^Fe levels out while ^^Fe actually decreases in abundance as the over- 
abundance of free neutrons drops. Meanwhile, as the free proton abundance rises toward 
equilibrium, ^^Ni and ^^Ni rise to prominence. This results in temporal variations in the 
relative abundances of the nuclei which dominate the distribution. It is not until very late 
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Figure 3a: Evolution of the mass fractions of tlie dominant members of the iron peak group as a function 
of the degree of siUcon exhaustion for Tg = 5.0, p = 10^ gcm""^, and = .498. Compare with the sohd line 
in Figure 2. 

that ^^Ni, which is the most abundant nucleus in the NSE distribution for these conditions, 
becomes most abundant in the network distribution. The relative abundances of ^^Fe, ^^Fe, 
and 5^Ni at X(Si group) ~ .9 are similar to those found in the NSE distribution for Yg ~ .475. 
At X(Si group) ~ .5, the relative abundances of ^^Fe, ^^Ni, and ^^Ni are comparable to those 
found in NSE for Yg ~ -49. In another similarity to the NSE distribution, nuclei with similar 
Z/A react in concert. For example the abundances of ^^Fe and ^*^Ni react to a decreasing 
silicon group mass fraction by rising sharply, peaking around X(Si group) ~ .9, and then 
falling off. For ^^Fe this shallow decline turns around for X(Si group) ~ .1, coincident with 
the minimum of the free neutron distribution. As the free neutrons converge back toward 
their equilibrium value the abundance of ^^Fe also increases slightly. For ^*^Ni this decline 
after maximum is turned around for X(Si group) ~ .4 as the increase in the free proton frac- 
tion overcomes the decline in the free neutron fraction. Even after the abundance of ^^Fe 
begins to increase again, the rate of increase in the abundance of ^°Ni, aided by its greater 
dependence on the free proton abundance, is greater. Other such pairings of nuclei observed 
in the NSE distributions, like ^^Fe and ^^Ni or ^^Fe and ^^Ni, also behave in concert. 

This excessive neutron-richness of the iron peak group at early times is also apparent 
for other thermodynamic conditions. For Tg = 5.0, p = lO^gcm""^, and Yg = -48, shown 
in Figure 3b (and as the dotted line in Figure 2), the relative distribution of ^^Fe, ^^Fe, 
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Figure 3b: Evolution of tlie mass fractions of tlie dominant members of the iron peak group as a function 
of the degree of sihcon exhaustion for Tg = 5.0, p = lO^gcm"^, and Yg = .48. Compare with the dotted 
Una in Figure 2. 



^'^Ni, and ^^Ni, at X(Si group) ~ .5, is more in keeping with equilibrium distributions with 
Yf> ~ .47. Once again, relatively neutron-rich isotopes, in this case ^^Fe, ^^Ni, and ^^Ni, peak 
early, exhibit a local minimum in concert with a similar local minimum in the abundance 
of free neutrons, and then rise toward equilibrium. Here again, the increasing abundance of 
free protons causes the Ni isotopes to rise more quickly than their Fe brethren. Thus the 
behavior of this example is consistent with the previous case, though even more neutron-rich, 
as one might expect from this lower case. For Tg = 3.5, p = lO^gcm^"^, and = -46, 
shown in Figure 3c (and as the dashed line in Figure 2), the early behavior is even more 
neutron-rich. For X(Si group) ~ .5, the relative distribution of ^^Fe, ^^Fe, ^°Ni, and ^^Ni 
reflects a Yg < -45. Even for material which is highly neutronized, the iron peak group is 
especially neutron-rich at early times. 

Over this wide range of conditions, the abundance distributions within the iron peak 
group exhibit a number of similarities. First, the relative abundances are very reflective 
of the changes in the free nucleon abundances. Second, at early times, particularly when 
the silicon group still dominates the mass fraction, the iron peak group exhibits a higher 
degree of neutronization than the abundance distribution as a whole. The solid line in Figure 
4, which plots the average Z/A of the iron peak group as the degree of silicon exhaustion 
increases, for Tg = 5.0, p = lO^gcm"^, and Yg = -48 (the same conditions as in Figure 
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Figure 3c: Evolution of the mass fractions of the dominant members of the iron peak group as a function 
of the degree of siUcon exhaustion for Tg = 3.5, p = lO^gcm"^, and = .46. Compare with the dashed 
Una in Figure 2. 

3b), shows that this excessive neutronization is not an illusion of the most abundant nuclei. 
For the iron peak group to show excessive neutronization there must be a group of nuclei 
which are neutron poor relative to the global Ye- This is required by the constancy of total 
neutron number in the absence of weak interactions. Since the mass fractions of the light 
particles are very small, the hkely source of these excess neutrons is the sihcon group. The 
dotted line in Figure 4, which shows the average Z/A for the silicon group as it is depleted, 
vindicates this assertion. Clearly, in equilibrium, the neutron-richness of the silicon group 
is much lower than the global neutronization indicates. While in small part this is due to 
the absence of the most neutron-rich isotopes of Ca, Sc, and Ti from the sihcon group, the 
major cause is the relative binding energies within the silicon group. 

As Table 3 indicates, under conditions of nuclear statistical equilibrium, even with — 
.48, elements which compose the silicon group are dominated by their N = Z isotopes. 
This is particularly true for the dominant elements. Si, S, Ar, and Ca, where the N = Z 
isotope is also an a-nucleus. As for nuclei with Z between 2 and 10, the valley of maximal 
binding energy is narrow for elements within the silicon group. With the high temperatures 
of silicon burning, and the resultant photon distribution, the more neutron-rich isotopes 
within the silicon group, with binding energies smaller than their N = Z neighbors, are 
more quickly photodissociated, freeing large numbers of free neutrons. These neutrons are 
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Figure 4: Evolution of the average Z/A for the iron peak group (solid line) and the silicon group (dotted 
line) with the degree of silicon exhaustion for Tg = 5.0, p = 10^ g cm~^, and Yg = .48. Compare with Figure 
3b. 



Table 3 Average Z/A function of Element 

in NSE at T., = 5.0. p - lO^gcm-^. and Y, - .18 



Element 


Z/A 


Element 


Z/A 


Element 


Z/A 


n 





Na 


.479 


Ti 


.477 


H 


1 


Mg 


.500 


V 


.473 


He 


.500 


Al 


.482 


Cr 


.474 


Li 


.500 


Si 


.500 


Mn 


.473 


Be 


.548 


P 


.485 


Fe 


.478 


B 


.464 
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.498 


Co 


.481 
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.500 


CI 


.486 


Ni 


.483 


N 


.500 


Ar 


.496 


Cu 


.480 





.500 


K 


.487 


Zn 


.480 


F 


.528 


Ca 


.499 


Ga 


.478 


No 


.500 


Sc 


.181 


Go 


.177 



captured by nuclei in the iron peak group, where the binding energy valley has curved to the 
neutron-rich side and widened, and thus the neutron-rich isotopes are more durable. Thus 
it is the tendency toward N = Z nuclei within the silicon group which fuels the extreme 
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neutron excesses in the iron peak group during silicon burning. This provides a mechanism 
for producing small quantities of more neutron-rich iron peak group nuclei, characteristic 
of a larger rj. As Figure 4 shows for Tg = 5.0, p = lO^gcm"^, and Yg = .48, if only 
10% of the silicon was exhausted, the effective rj within the iron peak group would be .08, 
twice its equilibrium value. Thus, by such a mechanism it is possible to produce neutron- 
rich members of the iron peak group from material which has never experienced the high 
densities necessary to produce such a large global neutron excess. By tying up the bulk of 
the nucleons in N=Z nuclei, and thereby amplifying the effective neutronization among the 
iron peak nuclei, incomplete silicon burning faintly echoes the emerging understanding of the 
r-process (see Woosley & Hoffman 1992). This neutron enrichment of the iron peak group is 
clearly a consideration when one is modeling objects where explosive burning of silicon, and 
hence incomplete silicon burning, are important to the elemental production. 

6. Reaction Flows 

With the majority of the mass concentrated within the two QSE groups, the process of 
silicon burning is dominated by these groups of nuclei. The previous section indicates that, 
with knowledge of the free nucleon abundances, it is possible to ignore the reactions within 
the quasi-equilibrium groups, since these reactions simply reflect the changes in abundances 
required for nuclei to remain in quasi-equilibrium under changes in the free-proton and free- 
neutron fractions. Thus, instead of studying the thousands of reactions which make up silicon 
burning, we can concentrate on those reactions which enter or leave the QSE groups. Ideally, 
these few important reactions would be the same for all physical conditions. However, in §5 
we demonstrated that small changes in neutronization result in much larger changes in the 
free-neutron abundance, thus favoring neutron capture and more neutron-rich nuclei. This 
enhanced importance of neutron captures affects the boundary region and results in changes 
in the location of the boundaries of the silicon and iron peak QSE groups, as we discussed in 
§4. As a result we expect that the dominant reactions linking these groups must change as a 
function of Yg. However, before we examine the complicated behavior that joins the silicon 
and iron peak groups, let us examine the lower boundary condition for the QSE groups. 

6.1. Downward Flows from Silicon 

While we have heretofore concentrated on the upper edge of the silicon quasi-equilibrium 
group, the flow from the lower edge is also important to the destruction of silicon and the 
other constituents of the silicon group. Previous authors have shown that these reactions 
govern the rate at which silicon is destroyed. Examination of Figure la, lb, or Id shows 
that the network abundances of the nuclei below A ~ 24 are orders of magnitude less than 
their silicon QSE abundances, with the degree of underabundance increasing as Z decreases. 
Further, there is little evidence of quasi-equilibrium behavior among these nuclei. Figure 
Ic and le reveal that the principal isotopes of Mg and Al are close to their silicon QSE 
abundances. Within a margin of ±10%, all of the isotopes of Mg and Al, except ^^Mg and 
perhaps ^^Mg, are members of the silicon group. This closeness to quasi-equilibrium does 
not hold for Na and Ne. Thus the bottom edge of the silicon group is Mg. If the principal 
isotopes of Mg and Al are in QSE with ^®Si, then the rate of photodissociation of ^^Si, 
either by ^®Si(7, a)^^Mg or by ^^Si(7, p)^^Al, is essentially balanced by the reverse capture 
reactions. Thus it is the net flow from this bottom edge which governs the downward flow 
from silicon. BCF contended that, with ^^Mg in quasi-equilibrium with ^^Si, the rate of 
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Figure 5a: Vector field representing the integrated reaction fluxes from X(Si group) =.95 to .5, along the 
lower boundary of the silicon group, for Tg = 5.0, p = lO^gcm"^, and Yf. = .498. The vector magnitudes 
have a logarithmic dependence on the size of the integrated flux. 

destruction of ^®Si was governed by the photodissociation of ^^Mg, with ^^Mg(7, a)^^Ne 
dominating ^^Mg(7, p)^'^Na. This was supported by WAC, who found that these reactions 
did indeed dominate the downward flow. However, the reaction network used by WAC 
was very narrow at this lower boundary of the silicon group with four isotopes of Mg, and 
only single isotopes of Al, Na, and Ne. Further, because of doubt over the accuracy of the 
then best available reaction rates, a number of reactions which those authors described as 
potentially important, like ^'*Mg(p, Q;)^^Na and ^'*Mg(n, a)'^^Nc, were excluded from their 
reaction network. Thus we now reexamine the reactions which cross this lower boundary to 
see the importance of our improved reaction network in addition to examining the effects of 
the range of Ye- 

While the variation of Yg has considerable influence on the relative abundances of the Mg 
isotopes, there is not the pronounced shift in group membership at this lower boundary that 
occurs at the upper boundary. Comparison of Figures Ic, If, and Ig, shows changes in the 
relative spacing of the Mg isotopes, and a somewhat more pronounced spread in those cases 
further from NSE, but little evidence of significant movement of the lower boundary. Thus 
the effects of variation of Yf, on the lower boundary are simply changes in relative abundances. 
Since changes in Fg have no counteracting effect on the reaction cross sections, these changes 
in relative abundances significantly alter the reaction balances. In this and the following 
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Figure 5b: Vector field representing the integrated reaction fluxes from X(Si group) =.95 to .5, along the 
lower bomidary of the silicon group, for Tg = 5.0, p = lO^gcm"^, and Yf. = .46. The vector magnitudes 
have a logarithmic dependence on the size of the integrated flux. 

subsection, we will show a series of vector diagrams of the integrated reaction fluxes, forward 

reaction rate minus its respective reverse, integrated over time. The direction of the vector 
reflects the dominant reaction of the pair, pointing from target to product. In order to show 
best the many reactions which contribute and yet not have the dominant vectors be visually 
overwhelming, the magnitudes of these vectors have a logarithmic dependence on the actual 
magnitude of the flux. As an example, in Figure 5a, for Tg = 5.0, p — lO^gcm"^ and Ye = 
.498, the flux through ^'*Mg(7, a)^^Ne is 12 times greater than that through ^■^Mg(n, a)^^Nc, 
although the vector representing ^^Mg(7, a)'^'^Nc is only 1.4 times longer than the vector 
representing ^^Mg(n, Q;)^^Ne. Figure 5a shows the reaction fluxes integrated from X(Si group) 
= .95 to .5, which is the portion of the burning where the boundary between the silicon group 
and elements of lower Z is most important. Even at X(Si group) = .5, the abundances of ^^Na 
and 20 Ne are less than ^ and ^ of their respective silicon quasi-equilibrium abundances. As 
expected, the downward flows from the silicon group strongly dominate mass transfer across 
this border. Only a few small fluxes, most notably 22Na(n, p)22Mg and 22Na(Q;, n)2^Al, 
point upward across the boundary. At first glance, this figure also seems to agree with the 
results of WAC, with the fiux through •^^Mg(7, a)20j\jg being nearly 3 times the fiux through 
2^Mg(7, p)2^Na. But the largest reaction fiux downward from the silicon group is actually 
2^Mg(p, Q;)2'^Na with a fiux 5.3 times larger than that through •^^Mg(7, p)2'^Na. In addition 
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Figure 5c: Vector field representing the integrated reaction fluxes from X(Si group) =.95 to .5, along the 
lower boundary of the silicon group, for Tg = 3.5, p = lO^gcm"^, and = .498. The vector magnitudes 
have a logarithmic dependence on the size of the integrated flux. 

there is a group of lesser flows, including ^^Mg(p, Q;)^^Na, ^■^Mg(n, p)^^Na, ^■^Mg(n, Q;)^^Ne, 
and ^^Al(n, Q;)^"^Na, which in sum rival the importance of ^'*Mg(7, a)^^Ne. This is certainly 
more complicated than the simple approach taken by BCF. Thus our more extensive network 
indicates that the apparent simplicity of the results of BCF and WAC is due to hmitations 
in their calculation. 

This is further reinforced by examination of the reaction fluxes for lower Yg. In Figure 5b, 
identical to Figure 5a, except that Yg = .46, the reactions ^^Mg(n, a)^^'Ne, ^^Mg(p, a)^"^Na, 
2^Mg(n,a)22Ne, 23Mg(n, p)23Na, '^^Mg{n,af^Ne, and 24Mg(n, p)24Na all carry more flux 
than 2^Mg(7, a)^^'Ne. This tendency at lower Y^ to favor reactions with neutrons in the 
incoming channel and neutron-rich nuclei as targets is very much in keeping with the much 
larger free-neutron abundances found for these conditions. There is also a noteworthy flux up- 
ward via 23Na(a, nf^Al This flux is actually part of a cycle, ^^Na{a, n)26Al(n, p)26Mg(p, af^Na§ 
with the flux through '^^Na.{a, n)^^A[ representing less than 7% of the flux through ^^Mg(p, a)2'^Na.| 
It is a myriad of cycles like this within the QSE groups which keep the a-particles in equi- 
librium with the free nucleons. 

Variations of temperature and density do not cause strong changes in the relative impor- 
tance of reactions like those due to variations of Yg. Figure 5c, with Tg = 3.5, p = 10^ g cm""^, 
and Yg = .498, is similar to Figure 5a. In this case 2^Mg(7, Q!)2^Ne does actually carry the 
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Figure 5d: Vector field representing the integrated reaction fluxes from X(Si group) =.1 to .02, along the 
lower boundary of the silicon group, for Tg = 5.0, p = lO^gcm"^, and Yf. = .498. The vector magnitudes 
have a logarithmic dependence on the size of the integrated flux. 

largest flux, but the same reactions we found important in Figure 5a are also important 
here. Thus the variation of temperature, which alters the balance of each set of forward and 
reverse reactions, and density, which enhances captures, can result in differences in relative 
fluxes, but in general does not open up drastically different paths, since these variations do 
not result in the large free neutron fraction found to be important for low Yg- With the 
passage of time and the approach to equilibrium, the reaction flows dwindle. The approach 
to equilibrium also implies that the isotopes of Na and Ne are not as under abundant. Thus 
the dominance of the flows downward from the silicon group over those flows directed upward 
lessens. This is illustrated in Figure 5d, which presents the integrated flux from X(Si group) 
= .1 to .02, for Tg = 5.0, p = lO^gcm"^ and Ye = -498. Over the corresponding interval 
of time, the abundance of ^^Ne converges from 80% of its QSE abundance to its NSE value. 
The sum of the upward fluxes, principally ^^Na(Q;, p)^^Mg, is almost as large as the sum 
of the downward fluxes. Furthermore, the balances between some pairs of reactions have 
changed. Representative of this trend, the balance of ^^Ne(Q;, n)^^Mg and ^^Mg(n, Q;)^-'^Ne, 
which favored ^^Mg(n, Q;)^-'^Ne at early times, favors ^^Ne(Q;, n)^^Mg as equilibrium is ap- 
proached. 

Thus the use of our larger and improved network indicates that there are a number of 
reactions, in addition to the photodisintegration reactions of ^^Mg, which are important to 
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understanding the depletion of the sihcon group toward Ughtcr nuclei. Under some con- 
ditions, particularly low Ye, these alternative reactions dominate the ^'^Mg(7, Q;)^^Ne and 
^^Mg(7, p)^^Na reactions which previous authors found to be the most important. We do 
agree with previous authors in finding that, as equilibrium is approached, the decreasing 
abundance of the Mg isotopes, mandated by quasi-equilibrium, and the closer proximity of 
the Ne and Na isotopes to their silicon quasi-equilibrium abundances, results in a decline 
in the net flux of mass downward from the silicon group. This equilibration of the lighter 
nuclei with the silicon group is the final phase of silicon burning, occurring well after the 
abundances are dominated by the iron peak elements and a single quasi-equilibrium group 
stretches upward from Mg. By this point, as we will show in §§6.2 and 7, the conversion of 
silicon into iron peak elements no longer dominates the energy production but serves as a 
source of free particles, driving the distribution to NSE. 

6.2. Bridging the Gap between Silicon and Iron 

While the downward fiow from the silicon group generates the free nucleons needed to 
build silicon into iron peak elements, and hence governs the timescale for silicon burning, 
the links between the silicon and iron peak quasi-equilibrium groups mediate the formation 
of the iron peak. For times prior to the establishment of a single QSE group, a small 
number of reactions dominate the fiow into the iron peak group. It is the slowness of these 
reactions which allows the persistence of two separate groups. WAC found that, with their 
network, the flow was dominated by a group of reactions ending in ^^Ti. They contended 
that the dominant flow was '*^Sc(p, 7)'*^Ti, aided at high temperature by ^^Ca(a, 7)'*^Ti 
and ^^Ti(n, 7)^*^Ti. This contradicted earlier work by BCF and Michaud & Fowler (1972), 
which had selected ^^Ti(Q;, p)^^V as the principal bridge. WAC argued that this reaction 
was only important at late times when a single QSE groTip was a good approximation. We 
will examine both of these contentions within the context of our larger network, and also 
examine the influence of Yg, something not done by previous authors. 

For Tg = 5.0, p = lO^gcm^"^, and Y^ = .498, we flnd that at early times there are 
two definite quasi-equilibrium groups separated by a thin boundary region (See Figure Ic). 
Although there are some relative changes in the abundances of these nuclei, variations of 
temperature, density, and degree of silicon exhaustion do not sharply alter the membership 
of this boundary region. In Figure 6a this boundary region is enclosed by the dotted line, 
with the lower edge of the iron peak group to the right and above, and the upper edge of the 
sihcon group to the left and below. As revealed in Table 2, for Tg = 5.0, p — 10^ g cm~^, and 
Ye = .498, by X(Si group) =.5 there is a single well established QSE group. In fact, for these 
conditions, the species which comprise the iron peak group have abundances within 10% 
of their silicon quasi-equilibrium abundances by the time X(Si group) ~ .85. As a result, 
the vectors in Figure 6a represent the integrated reaction fiuxes from X(Si group) =.95, by 
which time the two QSE groups are well estabhshed, until X(Si group) =.85. Thus these 
are the fluxes which bring the groups into quasi-equilibrium. It is therefore not surprising 
that there are no signiflcant fluxes downward from the iron peak group into the boundary 
region, and similarly no fiuxes into the silicon group which originate above. The dominant 
flux into the iron peak group is ^■'^Sc(p, 7)^^Ti. A number of other fluxes contribute including 
45Ti(n, 7)'^6rjj^ ^'2Ca{a,-f)'^^Ti, ^^Sc{a,p)'^^Ti, ^^V{n,p)^^Ti, and '^'^Ti{a,p)^^Y with fluxes 
relative to ^5Sc(p, 7)^6Ti of .12, .09, .04, .05, and .08, respectively. Thus '^5Sc(p, 7)^6Ti ac- 
counts for approximately 70% of the entire flux into the iron peak group. This agrees well 
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Figure 6a: Vector field representing the integrated reaction fluxes from X(Si group) =.95 to .85, along the 
boundary between the silicon group and the iron peak group, for Tg = 5.0, p = 10''gcm~^, and Yg = .498. 
The vector magnitudes have a logarithmic dependence on the size of the integrated flux. 

with the findings of WAC for explosive burning at comparable densities. This agreement 
occurs in spite of differences over group membership between this work and WAC. In par- 
ticular, WAC place ^^Sc as a member of the silicon group (those authors do not include a 
boundary region), while our results place it within the iron peak group for comparable Yg- 
The existence of the boundary region causes a further complication, in that there is not a 
one to one correspondence between the flows out of the silicon group and those into the iron 
peak group. In addition to reactions like '*'*Ti(«, p)^^V, '*"^Sc(a, p)^^Ti, and '^^Ca(a, 7)^^Ti, 
which directly connect the silicon group to the iron peak group, there are a number of reac- 
tions which transfer mass from the silicon group into the boundary region. Principal among 
these are '^'^Ca.{a,nf^Ti, 42Ca(Q;, p)45Sc, 44Ti(n, 7)4^X1, 44sc(p, 7)45x1, and ^^Scfn, p)44Ca, 
with fiuxes relative to 45Sc(p, 7)^6x1 of .28, .35, .21, .20, and .12, respectively. With this 
many reactions contributing significantly to the fiow of mass into the boundary region, the 
dominance of 45Sc(p, 7)4^X1 as a fiow out of the boundary requires substantive reactions 
among the nuclei of the boundary region. Thus the relative values of rqsEi'^Z) revealed 
in Figure Ic for the nuclei within the boundary are representative of reaction fiows within 
the region and not just a coincidental occurrence. Note in particular in Figure Ic that both 
^^Ti and ^^Ca are relatively overabundant when compared to "^^Sc and that Figure 6a shows 
significant fiow along 45xi(n, p)45Sc and 44Ca(p, 7)45Sc. Thus the links within the boundary 
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Figure 6b: Vector field representing the integrated reaction fluxes from X(Si group) =.75 to .1, along the 
boundary between the silicon group and the iron peak group, for Tg = 5.0, p = 10''gcm~^, and Yg = .498. 
The vector magnitudes have a logarithmic dependence on the size of the integrated flux. 



region arc important to the behavior of the boundary. The thin fringe of nuclei shown in 
Figure Ic between the sihcon and iron peak groups are not separate beads on strings of 
reactions between the QSE groups, but part of an interconnected region which mediates the 
transfer of mass between the sihcon and iron peak groups. 

Once these two QSE groups are replaced by a single group, knowledge of the flows 
between the formerly separate groups is no longer essential, as these flows represent the 
changes required under QSE as the free nucleon abundances converge toward equilibrium. 
In the case of Tg = 5.0, p = lO^gcm^"^, and Y,. = .498, this means that the fluxes for 
X(Si group) < .8 reflect the movement of the mass into more proton-rich nuclei, as we 
discussed in Sect. 5. Figure 6b shows the integrated fluxes from X(Si group) = .75 to .1. 
The dominant fluxes between the former silicon and iron peak groups, are a set of {a, p) 
reactions, principally '*'^Ti(Q;, p)^^V, ^^Sc(a, p)^^Ti, and '*^Ti(Q;, p)^^V. Thus the dominant 
mass fluxes at late time proceed upward through more proton-rich nuclei than they did prior 
to the merger of the QSE groups. Another reflection of this tendency toward proton-rich 
nuclei is the domination of (n,p) reactions by (p,n) after group merger. Comparison of 
Figures 6a and 6b reveals that, while at early time there was considerable flux within the 
iron peak QSE group through reactions like '^®Cr(n, p)^^V and '^^V(n, p)'*^Ti, after merger 
it is the reverse reactions, ^^V(p, n)'^^Cr and ^^Ti(p, n)^^V in this case, which dominate. In 
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Figure 6c: Vector field representing the integrated reaction fluxes from X(Si group) =.95 to .85, along the 
boundary between the silicon group and the iron peak group, for Tg = 3.5, p = lO^gcm"^, and Yg = .498. 
The vector magnitudes have a logarithmic dependence on the size of the integrated flux. 



another important reversal, at late times there are downward flows from the former iron peak 
group, most importantly '*^Ti(7, p)'*^Sc. For Tg = 3.5, p = lO^^gcm"^, and Y(, = .498, a 
considerably greater density than any considered by WAC, the scenario is little changed. As 
is shown in Figure 6c, the dominant fluxes, during the time that two separate QSE groups 
persist, are still 45Sc(p, 7)^6^1^ 45Ti(n, 7)46xi, 42ca(Q;, 7)46Ti, 43sc(o,, p)46Ti, 46v(p, n)46Ti, 
and ^^Ti(a, p)^'^V, with relative fluxes of 1, .05, .23, .007, .007, and .32, respectively. Clearly 
there is much change in the relative fluxes, with ^^Sc(p, 7)^^Ti carrying only 60% of the flux 
into the iron peak group, but the change is merely in degree. The variation of temperature 
and density does not profoundly change which reaction fluxes bridge the gap between the 
iron peak and silicon groups, although the relative importance can vary. Thus, for — .498 
our analysis of the dominant reaction flows which transfer mass into the iron peak group is 
in excellent agreement with that of WAC. The dominance of ^^Sc(p, 7)^^Ti as the primary 
conduit of mass transfer into the iron peak group is readily apparent. At late times the 
{a,p) reactions on more proton-rich nuclei dominate, especially ^^Ti(Q;, p)^^V, while there 
is a downward flow via '^^Ti(7, p)^^Sc. Both of these points support the earlier analysis of 
WAC. 

We noted earlier that the variation of Yg causes the membership of the QSE groups to 
change. These changes are reflected in the comparison of Figures Ic and le. As we discussed 
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Figure 6d: Vector field representing the integrated reaction fluxes from X(Si group) =.95 to .5, along the 
boundary between the silicon group and the iron peak group, for Tg = 5.0, p = 10''gcm~^, and Yg = -46. 
The vector magnitudes have a logarithmic dependence on the size of the integrated flux. 

in §4, the boundary region for smaller Ye is displaced and noticeably broader in N. This new 

boundary region is enclosed within the dotted lines of Figure 6d, which shows the integrated 
fluxes for Tg = 5.0, p = lO^gcm^"^, and Ye = .46, from X(Si group) =.95 to .5. Under 
these conditions the silicon and iron peak groups reach ±5% of merger for X(Si group) ~ .1. 
There is some ambiguity as to the group membership of the neutron-rich isotopes of K (^^K, 
^^K, ^^K, and ^^K), as they seem to be members of the silicon group for X(Si group) ~ .9 
with abundances approximately 60% of their silicon QSE abundance, while the members of 
the iron peak group are more than an order of magnitude underabundant (Figure le). For 
X(Si group) ~ .5 (Figure If), these K isotopes are still only approximately 80% of their 
silicon QSE abundance, while the iron peak group has converged to within a factor of 2. 
Examination of Figure 6d reveals that the reaction fluxes from these K isotopes into the 
neutron-rich isotopes of Ca are much smaller than the corresponding fluxes from Ca into Sc 
and are also much smaller than the flows into the boundary region from less neutron-rich 
isotopes of K and Ar. Thus their inclusion within the silicon group seems reliable, although 
imperfect. The dominant fluxes out of the silicon group are 38Ar(a,7)42Ca, 39Ar( a, n) Ca, 
40Ar(«,n)43ca, ^'^K{^,nf^Ga., and 43K(p, n)43Ca, with relative fluxes of .47, .34, 1, .1, and 
.30, respectively. There is also a small group of reactions, most importantly ^^Ar(Q;, n)^^Ca, 
42Ar(a,n)45cjj^ ^^K{a,n)^^^c, ^^K(p, n)^^^^^^ ^^K{^^nf^Ca., which directly link the 



33 




N (neutron number) 

Figure 6e: Vector field representing the integrated reaction fluxes from X(Si group) =.5 to .1, along the 
boundary between the silicon group and the iron peak group, for Tg = 5.0, p = lO^'gcm"^, and Y^. = .46. 
The vector magnitudes have a logarithmic dependence on the size of the integrated flux. 

silicon group to the iron peak group. However, the sum of these fluxes is less than 15% of 
the flux carried by ^^Ar(«, n)^"^Ca. These small fluxes are also only a minor contribution to 
the flux into the iron peak group, which is overwhelmingly dominated by ^^Ca(n, 7)^^Ca. 
The flux through this reaction is 2.3 times larger than the flux through '*°Ar(Q;, n)^^Ca and 
thus represents over 90% of the flux into the iron peak group. The flow within the iron peak 
group is dominated by a series of (n, 7) reactions between isotopes of Ca, out to "^^Ca. The 
flow upward from Ca is largely by (p,n) and (a, n) reactions. This transfer of mass into 
the iron peak nuclei via the neutron-rich isotopes of Ca is reminiscent of the exploratory 
results TA found for lower mass stars. In Figure 6c we examined the integrated fluxes for 
low temperature and high density, conditions one expects for core silicon burning in lower 
mass stars, but with Ye = .498, and found little difference between the low temperature/high 
density case and the high temperature/low density case, except for an increase in the relative 
importance of flows through the more proton-rich nuclei. Prom Figure 6d, we conclude that 
the flows through neutron-rich Ca which TA described are due to the greater neutronization 
found in the core of less massive stars. 

Examination of Figure 6e reveals that the trend noted earlier for less neutron-rich nuclei 
to increase in importance at later times is also true for Yg = -46. With merger of the QSE 
groups occurring for X(Si group) ~ .1 under these conditions, the integrated flux from X(Si 
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group) =.5 to .1 is still reflective of the flows which bring the groups to equilibrium. At these 
later times, the reaction ^^Ar(a, 7)^^Ca carries a larger flux than the reactions "^^Ar(a, n)'*^Ca 
and ^^Ar(Q;, n)^^Ca combined. Although the dominant flow into the iron peak group is still 
^^Ca(n, 7)^^Ca, comparison of Figure 6e to Figure 6d shows that the series of large fluxes 
through (n, 7) reactions on Ca isotopes is greatly diminished for ^^Ca and ^^Ca. Instead the 
flow proceeds upward at lower N. This process continues after group merger as the average 
Z/A within the iron peak group approaches Y^. Clearly the convergence of the average Z/A 
within the iron group toward Yg, as we discussed in §5, plays a role in making proton-rich 
nuclei more important as silicon is exhausted. With the iron group dominated at early 
times by nuclei with Z/A < Y^, the reactions through the neutron-rich isotopes are naturally 
enhanced by the increased abundance. As the abundance of free neutrons declines from the 
initial overabundance, the relative importance of the less neutron-rich nuclei, reflected by the 
average Z/A, grows. This is particularly true at low Yg, since the silicon group is dominated 
by nuclei with Z/A ~ .5. 

Thus the reactions which bridge the gap between the silicon group and the iron peak 
group reflect the underlying changes in abundance. The interconnection of these groups is 
also reflective of the group membership and of changes in the relative importance of these 
members. For high Y^, the dominant flow into the iron peak group prior to the merger of 
the groups is ^^Sc(p, 7)^^Ti. The flow which feeds ^^Sc exits the silicon group via a series 
of reactions involving ^^Ca, ^^Ca, "^"^Sc, ^"^Sc, ^^Ti and '^^Ti. For somewhat smaller Ye the 
flow into the iron peak group via ^^Sc(p, 7)^^Ti is still important, but it is now fed through 
reactions involving ^^Ca, ^"^Ca, and ^^Sc. As the degree of neutronization increases, the 
dominant flow between the groups naturally winds its way through more neutron-rich nuclei 
as the abundances of these nuclei become more important. For Yg = .46 this means that the 
merger is achieved via a flow that proceeds through a long series of (n, 7) captures on Ca. 
As time proceeds, as we showed in §5, the average Z/A within the iron peak group converges 
toward the global value of Yg- The network of reactions responds to/causes this greater 
abundance of less neutron-rich nuclei by increasing the importance of bridging reactions 
involving more proton-rich nuclei and by increasing the flux through (p,n) reactions within 
the iron peak group. It is the neutronization of material which dominates the behavior of 
the QSE groups, and therefore the reactions that merge the groups. Although the variation 
of Yg has httle effect on the cross sections and rate coefficients, the reaction flows across the 
boundary change in response to the changing abundances within the groups. Since these 
reactions, particularly those from the bottom of the silicon quasi-equilibrium group, govern 
the rate of conversion of silicon into iron peak nuclei, these differences in reaction flows will 
clearly be reflected in the rate of energy production. 

7. Energetics of Silicon Burning 

Before we examine the energy production from the actual network calculations, it is 
instructive to make some simple energetic arguments about silicon burning. If we take as the 
flnal distribution of each of our constant temperature, density, and Yg cases the corresponding 
screened NSE distribution and compare this with the initial distribution, we can calculate the 
amount of energy available from silicon burning. Table 4 shows the binding energies of the 
initial and NSE distributions and their difference for a few cases, all in ergs/g. Comparison 
of cases which differ only in Yg reveals that much more energy is available at lower Y^. 
This reflects the larger binding energy per nucleon of ^^Fe, ^^Fe, and ^^Ni in comparison to 
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Table 4 The Energy Reservoir 
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^^Ni, a fact which, for example, accounts for much of the light curve of supernovae. The 
comparison of cases which differ in temperature and density shows a greater availability of 
energy at lower temperatures and higher densities. The lower the temperature, the more the 
ex.p(B {"^ Z) / kT) term in the nuclear Saha equation dominates, causing a more prominent 
peak in the abundance distribution composed of the most bound nuclei. Similarly, a higher 
density also favors more massive nuclei which (up to Ni/Fe) have a larger binding energy. 
For example, for Ye = .498 and p ^ 10^°gcm~^, the change from Tg = 3.5 to Tg = 5.0 
causes a change in the mass fraction concentrated in isotopes of Mn, Fe, Co, and Ni, from 
.9996 to .996. For Ye = .498 and p = lO^gcm-S, the same change in temperature results 
in changes in this mass fraction from .9991 to .918. With as much as 6 times more energy 
available per silicon nucleus at low Ye, even if the rate of silicon depletion is the same in 
all cases, there can be tremendous differences in the energy generation and therefore in the 
hydrodynamic conditions which cause this process. This is particularly important as electron 
captures cause Ye to drop. 

The lack of a single dominating reaction in silicon burning makes prediction of the energy 
production complicated. In general, the rate of energy production per gram by a nuclear 
process, e is given in terms of the reaction rates, (reactions cm~^s~-^), by 

6=?^^^, (15) 
P 

where Qij is the Q- value for the reaction between i and j. In cases where a single reaction 
dominates the process the sum is unnecessary. Although it is possible, as was shown above. 
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Figure 7a: Rate of energy generation as a function of degree of silicon exhaustion, for Tg 
p = lO^gcm""^, with Yg varying from .498 to .46. 



5.0, and 



to calculate the energy store available for silicon burning, the plethora of possible governing 
reactions revealed in §6 makes estimation of the temporal behavior of this energy release 
difficult. For each reaction ij in silicon burning, there is a reverse reaction ki which is 
approaching equilibrium with it. Thus Equation (15) becomes 



{nj - rM{ij)) Q 



P 



(16) 



where Qki = —Qij- As the abundances approach equilibrium, e approaches zero, not nec- 
essarily because of the exhaustion of fuel but because each pair of reactions is balancing. 
While Table 4 shows the global value of S = J^^^ e dt, it says nothing about the interplay 
of the various r^-'s and rk£S. These, however, determine the temporal behavior of e. 

In §6 we concluded that there are significant differences in the relative importance of 
the many reaction flows which typify silicon burning. In particular, the variation of Yg 
has a profound effect on which reactions are the most important links between the quasi- 
equilibrium groups and which reactions are primarily responsible for flow downward from 
the silicon group. Examination of Figure 7a shows that these differences in reactions result 
in tremendous variation in e as a function of X(Si group) and Yg. These particular curves are 
for Tg = 5 and p = lO^gcm"^. While the behavior is similar in each case as equilibrium is 
approached, the initial stages, which are the most important from the point of view of energy 
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Figure 7b: Rate of energy generation as a function of degree of silicon exhaustion, for Tg = 3.5, and 
p = lO^^gcm"^, with Yf. varying from .498 to .46. 



production, vary greatly. Clearly the different reactions which dominate as Ye varies result in 
very different behavior. As we noted in §4, the lower Ye cases take less time to reach a similar 
degree of silicon exhaustion. Table 5 shows the elapsed time for several degrees of silicon 
exhaustion. Silicon burning at Ye = .498 is much slower than silicon burning at Ye = .46. 
With the elapsed time to X(Si group) ~ .9 almost an order of magnitude larger for high Yg 
and the energy available, the global £, much smaller, it is not surprising that e for 1^ = .498 
is more than an order of magnitude smaller than e for Yg — .46 at early times. Further, 
with the elapsed time between X(Si group) ~ .5 and .1 being more than 100 times larger for 
Ye = .498 than for Ye = .46, it is not surprising that e falls off more precipitously for higher Ye. 
These very different temporal evolutions in energy generation are then directly attributable 
to the very different fluxes seen in §6. Clearly simple approximations which multiply a 
temperature- and density-dependent formula for e by a correction dependent on Ye are ruled 
out. The behavior with changing Ye is dependent on the complex interaction of the different 
reactions out of the silicon group. Figure 7b shows that at high density/low temperature 
there is a similar tangle of curves as a function of Yg. Any successful approximation of 
the energy generation for silicon burning needs to take into account the complex behavior 
of the reactions out of the quasi-equilibrium groups. As we established in §§4 and 5, the 
reactions within the groups can be well understood by the use of quasi-equilibrium. The 
many reactions, shown in §6 to be important for the understanding of the flows among the 
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Table 5 Elapsed time verses Degree of Silicon Exhaustion 
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groups, are also the governing reactions here, particularly those reactions which transfer 
mass downward from the silicon group. The sheer number of important reactions, and their 
response to changes in Fg and degree of silicon exhaustion, preclude a simple analytic model. 
Instead, in Paper II, we will discuss a promising approximation to the network calculation 
of silicon burning employing quasi-equilibrium. For now, however, we will discuss a few 
additional general characteristics of the energy generation during silicon burning. 

One particularly noteworthy feature, common to all of the cases pictured in Figure 7a, 
is the tail off, with the energy generation actually becoming an energy depletion before the 
silicon group reaches its equilibrium value. This occurs once the iron peak nuclei dominate 
the mass fraction and thus the energy generation from the destruction of silicon group 
elements has dropped. The convergence of the free proton and free neutron fractions to their 
equilibrium values results in an increasing portion of the mass within the iron peak group 
being found in elements with higher Z and A. These nuclei are not as bound as the Fe 
and Ni isotopes which dominate the NSE distribution, hence their production is endoergic. 
Furthermore, the production of the free nucleons is also endoergic. As Yg decreases, the 
equilibrium mass fraction of these nuclei increases, requiring more energy to produce them 
and therefore causing the tail-off to begin further from equilibrium, i.e., while the energy 
generated from the conversion of silicon into iron peak nuclei is larger. This results in the 
increasing gap between the last positive value of the energy generation and the equilibrium 
silicon group mass fraction as Yg decreases. For Yg ~ .5, this effect is exacerbated by the 
transition as silicon exhaustion increases in which nucleus is most abundant. It is this 
transition, demonstrated in Figure 3a, from ^^Fe and 5'*Fe to the less bound ^^Ni which is 
responsible for the sharp dechne of e seen for larger values of Yg- 

Figure 8a shows the effect of temperature variations on e for p = 10^° g cm-^ and Yg — 
.46. The resemblance these curves share is striking. Each is a power law which tails off as 
each case approaches equilibrium. The lack of convergence between these curves indicates 
that changes due to variations in temperature are variations in rates, not changes in the 
relative importance of different reactions, which agrees with our analysis of the reaction 
fluxes in §6. From Table 4 it is apparent that the available energy supply is essentially 
unchanged by temperature variations for this density. Thus the differences in e are due to 
differing timescales for the burning. Comparison of the times necessary for these conditions 
to result in 10% exhaustion of the silicon group attest to this assertion. From the elapsed 
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Figure 8a: Rate of energy generation as a function of degree of silicon exhaustion, for p ■ 
Yg = .46, with Tg varying from 3.5 to 5.0. 
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times for X(Si group) ~ .9 shown in Table 5, we find that the differences in e are very 
similar to the differences in elapsed time. Clearly the increase in temperature increases the 
photodissociation rates, thus providing for quicker destruction of the silicon group. However, 
the variation in timescale with temperature is not consistent with the variations in the 
photodisintegration rate of ^^Si. This is in keeping with the argument by BCF, supported 
by our analysis, that it is not the photodissociation rate of silicon which governs the rate of 
silicon destruction; rather, it is the downward flow from the silicon group which determines 
the burning timescale. 

Variations in density have effects similar to those of temperature. Figure 8b portrays 
the effects of the variation of density, for Tg = 5.0, and Fg = -46. Once again we see power 
laws that tail off as each case approaches its equilibrium. Note that the variation of e with 
density, particularly at early times, is small. The variations in elapsed time to reach various 
degrees of silicon exhaustion as a function of density are also shown in Table 5. The ratio of 
e for p = lO^'^gcm""^ to that for p = lO^gcm""^ is only 1.5 for X(Si group) ~ .9, and 1.4 for 
X(Si group) ~ .5. Within the silicon QSE group, the increase in density results in an increase 
in the relative importance of Mg. The ratio of the mass fractions in Mg for p = 10^*^ gcm~^ 
compared to that for p = lO'^gcm"^ is 1.3 for both X(Si group) ~ .9 and X(Si group) ~ .5. 
This suggests that the difference in the energy generation rate as a function of density is 
largely due to differences in quasi-equilibrium abundances and in the reaction flows which 
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Figure 8b: Rate of energy generation as a function of degree of silicon exhaustion, for Tg = 5.0 and 
Yf, = .46, with p varying from 10''gcm~^ to lO^^gcm"^. 

result from them. 

In this section we showed that increases in temperature or density have large effects on the 
rate of energy generation, greatly enhancing the rate at which silicon is destroyed. Variations 
in neutronization also have a large effect on the rate of energy generation, changing the 
reaction pathways by which silicon is converted into iron peak elements. Thus, the evidence 
presented in this section points to the need to approximate the complete behavior of silicon 
burning in order to approximate the energy generation. This requires keeping track of many 
of the 300 nuclei and 3000 reactions used in this nuclear network. Fortunately, as we will 
show in a subsequent paper, quasi-equilibrium provides the means to approximate energy 
generation accurately at considerable savings in computation. Quasi-equilibrium greatly 
reduces the number of reactions which need to be considered, since only those which leave 
their group are important to the evolution of the abundances within the group. Further, 
the prediction of abundances based on the quasi-equilibrium abundances greatly reduces the 
amount of nuclear accounting which must be done. 

8. Conclusions 

We have performed a detailed reexamination of silicon burning as a function of temper- 
ature, density, and neutronization, using a large nuclear network. Central to understanding 
this process is the concept of quasi-equilibrium. A natural extension of nuclear statistical 
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equilibrium, which it approaches in the asymptotic hmit, quasi-equihbrium reflects the near 
balance of the rapid formation and destruction of nuclei during silicon burning. The net 
gain in abundance of a species is much smaller than either its formation or destruction rates. 
We rederived the equations of quasi-equilibrium, including for the first time the effects of 
Coulomb screening on the equilibrium distribution. This has been found by Hix et al. (1996) 
to be important for Nuclear Statistical Equilibrium, and we showed that Coulomb screening 
is also important in reconciling the network abundance calculations with quasi-equilibrium. 
We further demonstrated the real usefulness of quasi-equilibrium in describing the network 
evolution, with many of the most abundant nuclei forming two quasi-equilibrium groups, 
one focused on silicon and the other on the iron peak nuclei. While previous authors have 
seen such separate groups at early times, we discovered that, with increasing neutronization, 
the separation of these groups persists through a much more significant portion of silicon 
burning. Furthermore, we showed not only that the duration of this two group phase is a 
function of Y^, but also that the membership in the groups is dependent on Yf,. Since pre- 
vious work ignored extensive variation of neutronization and its effect on quasi-equilibrium, 
this behavior has not been seen before. 

Wc demonstrated the effects of these interacting quasi-equilibrium groups on the abun- 
dances as the distribution evolves. One noticeable effect is the excessive production of 
neutron-rich iron peak nuclei during incomplete silicon burning, resulting in an iron peak 
group noticeably more neutronized than the global neutronization would imply. This ef- 
fect is missed by many approximations to silicon burning, particularly those using narrow 
nuclear networks. With quasi-equilibrium governing the relative abundances within the 
groups, the principal evolution of the distribution is due to the reactions linking the groups. 
We further demonstrated that the important reactions for the destruction of silicon and the 
formation of iron peak nuclei are highly dependent on the degree of neutronization. While 
variations in temperature and density affect the relative rates of reactions and the overall 
speed of silicon destruction, they do not radically alter the reaction paths into and out of 
the quasi-equilibrium groups in the way that the variation of Yg does. For relatively unneu- 
tronized material, the principal reaction flow linking the silicon group with the iron peak 
group proceeds through isotopes of Sc and Ti on the proton-rich side of stability, in agree- 
ment with WAC. For larger neutronization, the increasing availability of free neutrons, and 
the corresponding changes in quasi-equilibrium group membership, result in the dominant 
flow proceeding through neutron-rich isotopes of Ca, as TA suggested. Furthermore, the 
rate at which matter flows along these reaction paths is also strongly enhanced under these 
conditions. With the entire process and especially the rate of silicon destruction strongly 
dependent on Y^, it is not surprising to find that the energy generation by silicon burning 
has a strong and complex dependence on Yg. As we also showed, the variations in energy 
generation due to changes in temperature and density are not as convoluted as those due to 
changes in Y^. 

Taken together, these effects of neutronization on the process of silicon burning imply 
that successful modeling of silicon burning in its hydrodynamic context needs to account for 
the variety of reaction paths and full range of important nuclei. Because large networks like 
the one used here are a cumbersome addition to hydrodynamic calculations, much previous 
work has tried to approximate silicon burning using narrow networks and other simplified 
schemes. The present results warn of the danger of applying such schemes, developed for 
slight neutronization {Y^, ~ .5), beyond the context in which they were developed. In the 
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context of reduced nuclear reaction networks, this implies that the networks not only need 
to stretch from H to the iron peak but also must have considerable breadth. This seemingly 
returns us to the cumbersome large networks. However, in a forthcoming paper we will 
demonstrate an alternative. Instead of reducing the number of nuclei by restricting the 
width, we find it possible to use the physics of quasi-equilibrium to intelligently reduce the 
number of abundances which must be evolved to a more manageable size. 

The hst of people who have contributed to this work is naturally larger than the author 
list. In particular, we would like to thank Dave Arnett, Ken Nomoto, and Masaki Hashimoto 
for motivating this work, and Al Cameron for his interest and aid. W.R.H. was supported 
in part by a NASA Graduate Student Researcher Fellowship. F.-K.T. was supported in part 
by NSF grant 89-13799 and the Swiss Nationalfonds. 
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